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1. Introduction 



One of the basic tasks of modem high energy physics is the study of the fundamental theory 
of strong interaction (Quantum Chromodynamics, QCD) in new, unexplored extreme regimes of 
super-high densities and temperatures through the investigation of the properties of multi-particle 
systems produced in high-energy nuclear collisions ^. Ongoing and future experimental heavy 
ion studies require the development of new Monte-Carlo (MC) event generators and improvement 
of existing ones. Especially for experiments at the CERN Large Hadron Collider (LHC), because 
of very high parton and hadron multiplicities, one needs fast (but realistic) MC tools for heavy ion 
event simulations ||3|,^ ||, ^. A realistic MC event generator should include a maximum possible 
number of observable physical effects which are important to determine the event topology: from 
the bulk properties of soft hadroproduction (domain of low transverse momenta pj < IGeV/c) 
such as collective flows, to hard multi-parton production in hot and dense QCD-matter, which 
reveals itself in the spectra of high-z^r particles and hadronic jets. However, in most of the avail- 
able MC heavy ion event generators, the simultaneous treatment of collective flow effects for soft 
hadroproduction and hard multi-parton in-medium production is absent. 

HYDJET-i-i- event generator ^] includes detailed treatment of soft hadroproduction as well 
as hard multi-parton production, and takes into account medium-induced parton rescattering and 
energy loss. The heavy ion event in HYDJET-i~i- is the superposition of two independent compo- 
nents: the soft, hydro-type state and the hard state resulting from multi-parton fragmentation. HY- 
DJET-i-i- model is the development and continuation of HYDJET event generator 10, 11, [I^ ], 



and it contains the important additional features for the soft component: resonance decays and more 
detailed treatment of thermal and chemical freeze-out hypersurfaces [|l3[ |l^]. The main program 
HYDJET-i-i- is written in the object-oriented C-i-i- language under the ROOT environment [15]. 



2. Physics model and simulation procedure 

The soft and hard components in HYDJET-i~i- are treated independently. When the generation 
of soft and hard components in each event at given b is completed, the event record (information 
about coordinates and momenta of primordial particles, decay products of unstable particles and 
stable particles) is formed as the junction of these two independent event outputs. 

The details on physics model and simulation procedure of HYDJET-i-i- can be found in the 
corresponding manual [0]. The main features of HYDJET-i~i- model are listed only very briefly in 
this section. 



2.1 Hard multi-jet production 

The model for the hard multi-parton part of HYDJET-i~i- event is the same as that for HYDJET 
event generator, and it based on PYQUEN partonic energy loss model [^, 1^, 11]. The approach 
to the description of multiple scattering of hard partons in the dense QCD-matter (such as quark- 
gluon plasma) is based on the accumulative energy loss via the gluon radiation being associated 
with each parton scattering in the expanding quark-gluon fluid and includes the interference effect 
(for the emission of gluons with a finite formation time) using the modified radiation spectrum 
dE/dl as a function of decreasing temperature T. The model takes into account radiative and 
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coUisional energy loss of hard partons in longitudinally expanding quark-gluon fluid, as well as 
realistic nuclear geometry. 



The Fortran routine for single hard nucleon-nucleon sub-collision PYQUEN |16] was con- 
structed as a modification of the jet event obtained with the generator of hadron-hadron interac- 
tions PYTHIA_6.4 [17]. The event-by-event simulation procedure in PYQUEN includes 1) gen- 
eration of initial parton spectra with PYTHIA and production vertexes at given impact parameter; 
2) rescattering-by-rescattering simulation of the parton path in a dense zone and its radiative and 
coUisional energy loss; 3) final hadronization according to the Lund string model for hard partons 
and in-medium emitted gluons. Then the PYQUEN multi-jets generated according to the binomial 
distribution are included in the hard part of the event. The mean number of jets produced in an 
AA event is the product of the number of binary NN sub-collisions at a given impact parameter 
and the integral cross section of the hard process in NN collisions with the minimum transverse 
momentum transfer p^f™. In order to take into account the effect of nuclear shadowing on parton 
distribution functions, the impact parameter dependent parameterization obtained in the framework 
of Glauber-Gribov theory is used. 

Note that some different approaches for MC treatment of partonic energy loss, such as codes 
YaJEM [|T9|], JEWEL ^ and Q-PYTHIA [|l]] have been presented during this Workshop. 



2.2 Soft "thermal" hadron production 

The soft part of HYDJET-i~i- event is the "thermal" hadronic state generated on the chemical 
and thermal freeze-out hypersurfaces obtained from the parameterization of relativistic hydrody- 
namics with preset freeze-out conditions (the adapted C-i-i- code FAST MC [ pj] , |l^). Hadron 
multiplicities are calculated using the effective thermal volume approximation and Poisson multi- 
plicity distribution around its mean value, which is supposed to be proportional to the number of 
participating nucleons at a given impact parameter of AA collision. The fast soft hadron simulation 
procedure includes 1 ) generation of the 4-momentum of a hadron in the rest frame of a liquid ele- 
ment in accordance with the equilibrium distribution function; 2) generation of the spatial position 
of a liquid element and its local 4-velocity in accordance with phase space and the character of mo- 
tion of the fluid; 3) the standard von Neumann rejection/acceptance procedure to account for the 
difference between the true and generated probabilities; 4) boost of the hadron 4-momentum in the 
center mass frame of the event; 5) the two- and three-body decays of resonances with branching ra- 
tios taken from the SHARE particle decay table [^2|]. The high generation speed in HYDJET-i~i- is 
achieved due to almost 100% generation efficiency of the "soft" part because of the nearly uniform 
residual invariant weights which appear in the freeze-out momentum and coordinate simulation. 

Let us indicate some physical restrictions of the model. HYDJET-i~i- is only applicable for 
symmetric AA collisions of heavy (A > 40) ions at high energies (y^ > 10 GeV). Since the hydro- 
type approximation for heavy ion collisions is considered to be valid for central and semi-central 
collisions, the results obtained for very peripheral collisions (with impact parameter of the order 
of two nucleus radii, b ~ 2Ra) may be not adequate. Nor do we expect a correct event description 
in the region of very forward rapidities, where the other mechanisms of particle production, apart 
from hydro-flow and jets, may be important. 
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particles.data, 
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Figure 1: The block structure of HYDJET++. 



3. HYDJET++ software structure 

The basic frameworks of HYDJET++ are preset by the object-oriented C++ language and 
the ROOT environment [15]. There is also the Fortran- written part Jl2[ ] which is included in the 
generator structure as a separate directory. The block structure of HYDJET++ is shown in Figure 
The main program elements are particle data files, input and output files, C++ and Fortran routines. 

The information regarding the particle species included in the HYDJET++ event is stored in 
the files particles . data and tabledecay . txt. The particles . data file contains the 
definition (PDG code) and physical properties (mass, decay width, spin, isospin, valence quark 
composition) of 360 stable hadrons and resonances. The tabledecay .txt file contains decay 
channels and branching ratios. The structure of these files is the same as that in SHARE particle 



data table [g2p and in event generator THERMINATOR [^]. 

Run of HYDJET++ is controlled by the file Run Input Hyd jet for different type of input pa- 
rameters. The input file contains 7 input parameters (number of events to generate, beam c.m.s. en- 
ergy per nucleon pair in GeV, atomic weight of nuclei and parameters to specify the type of central- 
ity selection) and 18 free model parameters, which can be varied by the user from their default val- 
ues (chemical potentials, chemical and thermal freeze-out temperatures, space-time scales of soft 
hadron emission, maximal longitudinal and transverse flow rapidities, momentum and coordinate 
azimuthal anisotropy parameters, PYQUEN energy loss model parameters). A number of important 
PYTHIA parameters also may be changed/specified in RunlnputHyd jet file. There are also 10 
flags to specify different physical model scenarios for soft and/or hard components. In particular, 
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the activation of some flags allows one to calculate a few model parameters, and so to reduce the 
number of independent free parameters (down to 12 at the minimum). Two additional files with the 
optimized parameters for Au+Au collisions at ^/s = 200A GeV (RunlnputHyd jetRHIC200) 
and for Pb+Pb coUisions at ^/s = 5500A GeV (RunlnputHyd jetLHC550 0) ai-e available in the 
distribution package. The default parameters for RunlnputHyd jetRHIC200 were obtained by 
fitting RHIC data to various physical observables. The default parameters for 
RunlnputHyd jetLHC5500 represent our rough extrapolation from RHIC to LHC energy. 

The program output is directed to the ROOT file RunOutput . root. The output file con- 
tains a tree named td, which keeps the entire event record including primary particles and decay 
products with their coordinates and momenta information. Each decay product contains the unique 
index of its parent particle so that the entire event history may be obtained. Beside particle informa- 
tion, the output file contains also a number of global output parameters for each event (generated 
value of impact parameter, total inelastic and hard scattering cross sections, hadron event multiplic- 
ities of hard and soft components, numbers of binary NN sub-collisions and nucleons-participants). 

HYDJET-i~i- includes 17 C-n- source files (4 main modules and 13 service modules placed in 
the main directory) and 3 Fortran files (placed in the separate directory). The size of the distribution 
package is 3.5 MBytes directory and 800 kBytes compressed archive (without ROOT libraries). 
The generation of 100 central (0-5%) Au-i-Au events at ^ = 200A GeV (Pb-i-Pb events at ^ = 
5500A GeV) with default input parameters takes about 7 (85) minutes on a PC 64 bit Intel Core 
Duo CPU @ 3 GHz with 8 GB of RAM memory under Red Hat Enterprise. Then the output file 
created by the code in ROOT tree format will require 40 (190) MBytes of the disk space. 



4. Validation of HYDJET++ with experimental RHIC data 



It was demonstrated in \ \13[ |14| ] that FAST MC model can describe well the bulk properties of 
hadronic state created in Au-i-Au collisions at RHIC at ^/s = 200A GeV (such as particle number 
ratios, low-pr spectra, elliptic flow coefficients V2{pT,b), femtoscopic correlations in central colli- 
sions), while HYDJET model is capable of reproducing the main features of jet quenching pattern 
at RHIC {high-pT hadron spectra and the suppression of azimuthal back-to-back correlations) |^]. 
Since soft and hard hadronic states in HYDJET-i~i- are simulated independently, a good description 
of hadroproduction at RHIC in a wide kinematic range can be achieved, moreover a number of 
improvements in FAST MC and HYDJET have been done as compared to earlier versions. A num- 
ber of input parameters of the model can be fixed from fitting the RHIC data to various physical 
observables. 



1. Ratio of hadron abundances. It is well known that the particle abundances in heavy ion 
collisions in a wide energy range can be reasonable well described within statistical models 
based on the assumption that the produced hadronic matter reaches thermal and chemical 
equilibrium. The thermodynamical potentials fls = 0.0285 GeV, /I5 = 0.007 GeV, Hq = 
—0.001, the strangeness suppression factor )J = 1, and the chemical freeze-out temperature 
T^^ = 0. 165 GeV have been fixed in [ 13] from fitting the RHIC data to various particle ratios 
near mid-rapidity in central Au-i-Au collisions at ^/s = 200A GeV (n^ /n^, p/n^ , /K^ , 
K-/TI-, p/p, A/A, A/A, S/S, ^/K-, K/p, E"/?!-). 
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2. Low-jcr hadron spectra. Transverse momentum pr and transverse mass my hadron spectra 
(71+, K+ and p with mt < 0.7 GeV/c^) near mid-rapidity at different centralities of Au+Au 



collisions at ^/s = 200A GeV were analyzed in ^1%. The slopes of these spectra allow the 
thermal freeze-out temperature = 0.1 GeV and the maximal transverse flow rapidity in 
central collisions p™''(^ = 0) = 1.1 to be fixed. 

3. Femtoscopic correlations. Because of the effects of quantum statistics and final state inter- 
actions, the momentum (HBT) correlation functions of two or more particles at small relative 
momenta in their c.m.s. are sensitive to the space-time characteristics of the production pro- 
cess on the level of fm. The space-time parameters of thermal freeze-out region in central 



Au-i-Au collisions at ^/s = 200A GeV have been fixed in [ 14 ] by means of fitting the three- 
dimensional correlation functions measured for 71+71+ pairs and extracting the correlation 
radii /?side, Rout and /?iong: Tfib = 0) = 8 fm/c, ATf{b = 0) = 2 fm/c, Rf{b = 0) = 10 fm. 



4. Pseudorapidity hadron spectra. The PHOBOS data on T] -spectra of charged hadrons [ |24[ ] 
at different centralities of Au-i-Au collisions at = 200A GeV have been analyzed to fix 
the particle densities in the mid-rapidity region and the maximum longitudinal flow rapidity 
'?max = 3.3 (Fig. HI). Since mean "soft" and "hard" hadron multiplicities depend on the cen- 
trality in different ways (they are roughly proportional to Npart{b) and Nbinib) respectively), 
the relative contribution of soft and hard parts to the total event multiplicity can be fixed 
through the centrality dependence of dN/drj. The corresponding contributions from hydro- 
and jet-parts are determined by the input parameters H^^^ = 0.053 GeV and p™" = 3.4 
GeV/c respectively. 

5. High-pT hadron spectra. High transverse momentum hadron spectra (pj > 2 — 4 GeV/c) 
are sensitive to parton production and jet quenching effects. Thus fitting the measured high- 
Pt tail allows the extraction of PYQUEN energy loss model parameters. We assume the 
QGP formation time To = 0.4 fm/c and the number of active quark flavours Nf = 2. Then 
the reasonable fit of STAR data on high-z^r spectra of charged pions at different centralities 



of Au-i-Au collisions at = 200A GeV [ 25 ] is obtained with the initial QGP temperature 
Tq = 0.3 GeV (Fig. |). 

6. Elliptic flow. The elliptic flow coefficient V2 (which is determined as the second-order 
Fourier coefficient in the hadron distribution over the azimuthal angle (p relative to the reac- 
tion plane angle that V2 = (cos2((p — Yr))) an important signature of the physics 
dynamics at early stages of non-central heavy ion collisions. According to the typical hy- 
drodynamic scenario, the values V2{pt) at low-px (< 2 GeV/c) are determined mainly by 
the internal pressure gradients of an expanding fireball during the initial high density phase 
of the reaction (and it is sensitive to the momentum and azimuthal anisotropy parameters 
5 and e in the frameworks of HYDJET-i-i-), while elliptic flow at high-;??- is generated in 
HYDJET-i-i- (as well as in other jet quenching models) due to the partonic energy loss in an 
azimuthally asymmetric volume of QGP. Figure ^ shows the measured by the STAR Col- 
laboration transverse momentum dependence of the elliptic flow coefficient V2 of charged 
hadrons in Au-i-Au collisions at = 200A GeV for two centrality sets p3]. The values of 
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Figure 2: The pseudorapidity distribution of charged hadrons in Au+Au collisions at ^/s = 200A GeV for 
six centrality sets. The points are PHOBOS data [|4|, histograms are the HYDJET++ calculations (solid - 
total, dotted - hydro part, dashed - jet part). 




[GeV/c] 



Figure 3: The transverse momentum distribution of positively charged pions in Au+Au collisions at 
= 200A GeV for three centrality sets. The points are STAR data histograms are the HYDJET++ 
calculations (solid - total, dotted - hydro part, dashed - jet part). 



5 and e are estimated for each centrality. Note that the choice of these parameters does not 
affect any azimuthally integrated physics observables (such as hadron multiplicities, rj- and 
pr-spectra, etc.), but only their differential azimuthal dependences. 
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Figure 4: The transverse momentum dependence of the elHptic flow coefficient V2 of charged hadrons in 



Au+Au colHsions at ^/s = 200A GeV for two centraHty sets. The points are STAR data |26|, histograms are 
the HYDJET++ calculations. 

5. Simulations with HYDJET++ at LHC 

The heavy ion collision energy at LHC a factor of 30 larger then that at RHIC, thereby al- 
lows one to probe new frontiers of super-high temperature and (almost) net-baryon free QCD. The 
emphasis of the LHC heavy ion data analysis (at ^/s = 5.5 TeV per nucleon pair for lead beams) 
will be on the perturbative, or hard probes of the QGP (quarkonia, jets, photons, high-pr hadrons) 
as well as on the global event properties, or soft probes (collective radial and elliptic flow effects, 
hadron multiplicity, transverse energy densities and femtoscopic momentum correlations). It is ex- 
pected that at LHC energies the role of hard and semi-hard particle production will be significant 
even for the bulk properties of created matter. HYDJET-i-i- seems to be an effective simulation tool 
to analyze the influence of in-medium jet fragmentation on various physical observables. 

Figures || and |6| show the pseudorapidity distribution of charged hadrons and transverse mo- 
mentum distribution of pions respectively obtained with HYDJET-i~i- default settings (in particular, 
p™" = 7 GeV/c) for 5% most central Pb-i-Pb events. The estimated contribution of hard component 
to the total event multiplicity is on the level ~ 55% here, what is much larger than as compared 
with RHIC (~ 15%, see Fig. §). Of course, this number is very sensitive to the parameter /j^'" 
— minimal pj of "non-thermalized" parton-parton hard scatterings. For example, increasing the 
value p™" up to 10 GeV/c results in decreasing this contribution down to 25%. 

Some applications of HYDJET-i~i- for high-pr studies at the LHC have been presented dur- 



ing this Workshop [E7L 28 1. In this paper we discuss one another striking example: the influence 



of jets on femtoscopic momentum correlations (HBT-radii) [29]. Since HYDJET-i-i- specifies the 
space-time structure of a hadron emission source (for soft and for hard components as well), the 
momentum correlation function can be introduced by the special weighting procedure [30, pT| ]. 



Knowing the information on final particle four-momenta and four-coordinates Xi of the emis- 
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Figure 5: The pseudorapidity distribution of charged hadrons in 5% most central Pb+Pb collisions at \/s ■ 
5500A GeV (soUd - total, dotted - hydro part, dashed - jet part), p??^" = 7 GeV/c. 
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Figure 6: The transverse momentum distribution of pions in 5% most central Pb+Pb collisions at \/s ■ 
5500A GeV (soUd - total, dotted - hydro part, dashed - jet part), pj?" = 7 GeV/c. 
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Figure 7: The ;r*;r* correlation radii and the strength parameter X at mid-rapidity as the function of relative 
pion momentum k in 5% most central Pb+Pb collisions at ~ 5500A GeV (black triangles - p™" — 7 
GeV/c, red triangles - = 10 GeV/c, blue triangles - pure hydro). The open circles are STAR data [^. 

sion points allows one to calculate the correlation function with the help of the weight procedure, 
assigning a weight to a given particle combination accounting for the effects of quantum statistics: 

w = 1 +cos(^- Ajc), (5.1) 

where q = pi — pi and Ax = xi — X2. Then the correlation function is defined as a ratio of the 
weighted histogram of the pair kinematic variables to the unweighted one. The corresponding 
correlation widths are parameterized in terms of the Gaussian correlation radii /?,, 

CF{pi,p2) = 1 + Aexp(-/?Qm^Qm -/?^ijig^^jde --^Lg^Lg -2/?Qm_iong^out^iong) , (5.2) 

where q = (^outi^side^^iong) is the relative three-momentum vector of two identical particles and A 
is the correlation strength. The out and side denote the transverse, with respect to the reaction axis, 
components of the vector q; the out direction is parallel to the transverse component of the pair 
three-momentum. 

The correlation functions of two identical charged pions have been calculated for 5% most 
central Pb-i-Pb events at two values of 7 GeV/c and 10 GeV/c. Figure ^ shows the corre- 
sponding values of fitted correlation radii and strength parameter as a function of relative pion 



momentum in comparison with those measured by STAR collaboration [32]. One can see that the 



"pure hydro" scenario results in some increasing of correlation radii at LHC in comparison with 
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RHIC, as it could be naively expected from the appropriate extrapolation of the volume param- 
eters However increasing the contribution of hard component to the total event multiplicity 
results in decreasing the correlation radii due to the fact that "jet-induced" hadrons are emitted on 
much shorter space-time scales than "thermal" particles. HYDJET-i-i- predicts that the correlation 
radii at LHC become comparable with those at RHIC for 25% contribution of hard component 
(p™'" = 10 GeV/c), and may be even less than at RHIC for larger contribution of hard component 
(p™'" < 10 GeV/c). On the other hand, the influence of hard component on correlation radii at 
RHIC was found to be negligible. Thus the observation of reducing the correlation radii as one 
moves from RHIC to LHC could manifest the strong influence of in-medium jet fragmentation on 
the bulk event topology. 

6. Summary 

Among other heavy ion event generators, HYDJET-i-H- focusses on the detailed simulation of 
jet quenching effect basing on the partonic energy loss model PYQUEN, and also reproducing the 
main features of nuclear collective dynamics by fast (but realistic) way. The final hadron state in 
HYDJET-i~i- represents the superposition of two independent components: hard multi-parton frag- 
mentation and soft hydro-type part. The main program is written in the object-oriented C-i-i- lan- 
guage under the ROOT environment. This model is the development and continuation of HYDJET 
event generator. The hard part of HYDJET-i-i- is identical to the hard part of Fortran- written HYD- 
JET and it is included in the generator structure as a separate directory. The soft part of HYDJET-i~i- 
is the "thermal" hadronic state generated on the chemical and thermal freeze-out hypersurfaces ob- 
tained from the parameterization of relativistic hydrodynamics with preset freeze-out conditions. It 
contains the important additional features as compared with HYDJET: resonance decays and more 
detailed treatment of thermal and chemical freeze-out hypersurfaces. HYDJET-i-i- is capable of 
reproducing the bulk properties of heavy ion collisions at RHIC (hadron spectra and ratios, radial 
and elliptic flow, femtoscopic momentum correlations), as well as high-pj- hadron spectra. 

HYDJET-i~i- is an effective simulation tool to analyze the influence of in-medium jet frag- 
mentation on various physical observables at the LHC. In particular, the spectacular prediction of 
HYDJET-i~i- is reducing the femtoscopic correlation radii in heavy ion collisions as one moves 
from RHIC to LHC due to the significant contribution of (semi-)hard component to the space-time 
structure of the hadron emission source. 
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